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Abstract: Recently we proposed that the original Berkovits-Witten twistor-string theory 
could be used to compute Einstein supergravity tree amplitudes. Subsequently, remarkable 
formulae for the MHV amplitudes were found by Hodges and generalized to N'^MHV by 
Cachazo and Skinner. Here we show how the Hodges matrix and its higher MHV-dcgrce 
generalization can be understood as arising from the twistor-string formula. In particular 
we see that if all contractions in the worldsheet correlator calculations are allowed, the 
wrong answers are obtained. On the other hand we obtain the correct answer if the con- 
tractions form connected trees (alternatively, the contractions in the correlators must leave 
a minimal number of disconnected terms) . The generalised Hodges matrix arise as weighted 
Laplacian matrices for the graph of possible contractions. The reduced determinants of 
these weighted Laplacian matrices give the sum of of the connected tree contributions by 
the Matrix-Tree theorem. 
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1 Introduction 

The twistor-string introduced by Witten [1] and realized as a momentum space tree formula 
[2] for AT = 4 super- Yang-Mills has been a remarkable stimulus for developments in our 
understanding of gauge theory amphtudes (see [3] for a review biased towards the interests 
of these authors and [4] and other reviews in that volume for other influences) . The question 
naturally arises as to whether analogous ideas can be made to work for gravity. Witten's 
twistor-string theory does contain conformal gravity, a fourth order conformally invariant 
theory whose action is given by the square of the Weyl tensor. However, at this point 
no twistor-string theory is known that produces just Einstein gravity, although the recent 
formulae of Cachazo and Skinner [5, 6] (see also [7, 8]) for J\f = 8 supergravity amplitudes 
from rational curves in twistor space are remarkably suggestive of the existence of such a 
theory. 

In [9] it was argued, using an observation of Maldacena [10], that the twistor-string 
tree formula for conformal supergravity amplitudes due to Berkovits and Witten [11] could 
be used to calculate Einstein tree amplitudes. The procedure yields an extra factor of the 
cosmological constant A and hence should vanish as A — t- 0. However, for an n-particle 
amplitude, the procedure automatically yields a polynomial of degree-n in A so that it is 
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straightforward to divide by A and take the Hmit A — t- 0. In [9], only three-point MHV 
and anti-MHV amphtudes were checked in momentum space (found for ah A and checked 
in the A — > hmit). In this article we investigate more general formulae that are yielded 
in momentum space in the A = limit. 

The first thing we observe is that the formulae in [9] must be understood in a semi- 
classical context in which only connected trees^ arc allowed amongst the contractions in the 
correlator in order to obtain the correct answer when A = 0. Indeed, the answer is simply 
non- vanishing at A = if loops and more disconnected terms amongst the contractions in 
the correlation functions are allowed. 

We move on to consider the 0(A) part of the conformal gravity amplitude, which by 
the Maldacena argument above yields the Einstein gravity contribution at A = 0. Working 
with just connected trees, we see that the tree formulae for the MHV amplitude of [12] can 
be identified simply with the Feynman diagrams of contractions required for the worldsheet 
correlators restricted to being connected trees. Furthermore, by summing the diagrams 
using a weighted extension of the Matrix-Tree theorem, we can obtain a more fundamental 
understanding of the origin of Hodges' recent remarkable MHV formula [13], as well as one 
of the reduced determinant factors in the N'^MHV formula of Cachazo and Skinner [5, 6]. 
These Hodges matrices can be understood as weighted Laplacian matrices for all possible 
contractions extended to a permutation invariant framework. The reduced determinants 
that arise in these formulae come from extensions of the Matrix- Tree theorem. We therefore 
obtain not only a more fundamental explanation for the origin of the ingredients in Hodges' 
formulae, but also support for the assertions of [9] for MHV, at least at 0(A) and with 
the proviso on the restrictions of allowable contractions in the twistor-string correlation 
functions. 

The conformal gravity twistor-string formulae are very much rooted in A/" = 4 super- 
symmetry and so do not manifest full permutation symmetry. However, if they are correct, 
they must have an emergent permutation symmetry presumably arising from a A/" = 8 
formulation of the momentum space analogue of the tree-formula. As already mentioned, 
such a formula has independently been found by Cachazo and Skinner [5, 6]. In this paper 
we see how the Hodges formula and certain key ingredients of the Cachazo- Skinner for- 
mula naturally arise from the Berkovits-Witten twistor-string theory. We obtain a fairly 
complete picture at MHV but rather less complete at higher MHV-degree. 

We now give a brief review of the Cachazo-Skinner formula. It is a natural extension 
of Hodges' formula for the MHV amplitude, and is a function of the kinematic invariants of 
the momenta AjAj where i = 1, . . . ,n index the external particles and Aj is a two-component 
Weyl spinor Ajyi, A = 0,1. For Aj, we understand it to also include the fermionic momenta 
so that Aj = (Ajyi', ?7") where A' = O'l' are spinor indices and o = 1, ... 8 is a i?-symmetry 
index. It also uses auxilliary variables CTj = {aio, an) that are homogeneous coordinates for 
n marked points on the Riemann sphere (the string worldsheet). We use square brackets 
for primed spinor contractions, angle brackets for unprimed spinor contractions and round 

^Rather than ruling out loops, we can require simply that the number of disconnected components 
amongst the contractions is minimal; we thank Mat Bullimore for this observation. 
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brackets for contractions of the as. 

Following [5, 6, 13], we introduce a n x n matrix function of the kinematic invariants 
and (TjS at N'^MHV de gree (i?-charge k + 2) by 



6^- 



fc+2 , (1-1) 



where ^ is an arbitrary fixed point on CP^. The off-diagonal entries of (f)'' are just the 
twistor-string propagators associated to contractions among worldsheet vertex operators. 
This is related to the matrix <I>'^' appearing in [5, 6] (and in [13] for k = 0) via conjugation 
by the matrix T = diag ((^ 



I P\ i^j 

^ = { ^ .;i/r.nX^+2 . 4>' = T-'^'T. (1.2) 



Hence, reduced determinants of (j)^ and will be equivalent^. 

When k = 0, (1.2) is Hodges' matrix for the MHV formula, which has three-dimensional 
kernel as a consequence of momentum conservation. The general has a A;-|-3-dimensional 
kernel given by the relations 

n 

which follow from the support of delta functions in the formula. A type of Fadcev-Popov ar- 
gument leads to an invariant determinant det'(^'^) of this matrix generalizing that described 
by Hodges. This is the main new ingredient in the tree formula beyond the Yang-Mills 
case. The final formula can be expressed as 

^nih...,n) = j ^^[fi72^det'($^)det'($'=) J] ^<^i hi (a,A,; Z{a,)) , (1.4) 

where the Ur are parameters for the degree A;-|-l map Z : CF^ PT; det'($'^) and det'($'^) 
are modified determinants; and the hi are momentum eigenstates. 

It is instructive to compare this to the analogous formula for A/" = 4 Yang-Mills: 

As remarked in [2], after integrating out part of the moduli space, both formulae have 
as many delta functions as integrals and the result reduces to a sum of residues at the 
solutions (both real and complex) of the equations implied by delta-functions. The main 



^The fact that (1.2) and hence (1.1) is independent of the choice ^ e CP^ follows from the delta^function 
support of the amplitude; see [5, 6] for more details. 
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difference between (1.4) and (1.5) is that the Parke- Taylor-like denominator of A/^ = 4 is 
replaced by the product of determinants det'(^>'^') det'($'^). 

This paper gives an interpretation of the det'(i>'^) factor in (1.4) from twistor-string 
theory. Indeed, the individual off-diagonal entries of the matrix are precisely all the possible 
propagators arising from a certain type of contraction in the world-sheet conformal field 
theory extended to a permutation invariant framework. The diagonal entries are then 
precisely what is required to obtain the weighted Laplacian matrix for the graph of all 
possible contractions. We will see that at N'^MHV, only n — 3 — k such contractions should 
be allowed and the reduced determinant is constructed precisely as the n — 3 — k minor 
that effectively sums trees with just n — 3 — k such propagators. 

There arc nevertheless significant gaps in our understanding here. In particular, we do 
not yet understand how the twistor-string ingredients can be used to assemble det'$, nor 
the Vandermonde factors in the definition of det' <I>. If the Berkovits-Witten twistor-string 
theory is to correctly describes tree-level conformal gravity, it must also be possible to 
obtain these within that framework. Indeed it is so closely connected to the twistor-action 
for conformal gravity [14], which has been systematically shown to give rise directly to 
classical conformal gravity, that it seems unlikely it will not be possible to obtain these 
ingredients, but a deeper understanding is required. Of course, of independent interest, 
the considerations in this paper arc also likely to play a role in any yet to be understood 
Af = 8 twistor-string theory for Einstein gravity. 

2 The Twistor-string Formula for Conformal and Einstein Gravity 

Here we briefly summarize the ingredients used in our twistor-string construction; more 
details can be found in [9]. Note that we abuse notation slightly, since in this section 
there is a SU(4)ij i2-symmetry associated with the twistor-string construction. We will use 
similar notation for both this J\f = 4 and the M = 8 i?-symmetry used elsewhere. 

2.1 Twistor-strings for conformal supergravity 

Non-projective twistor space is T = C^l^ and projective twistor space is PT = CP^'^ = 
T/{Z ~ e°^}, for a G C. A twistor wih be represented as eT, = {Z°',x"'), « = 
0, . . . 3, a = 1, . . . , 4 with .Z" bosonic and x"" fermionic and the bosonic part Z"' = {\a, ), 
A = 0,1, A' = O'l'. A point {x,9) = (x"^^' ,9^°-) in chiral super Minkowski space-time M 
corresponds to the CP^ (complex line) X C PT via the incidence relation 

fi^' =ix^'''XA, x'' = 0^''Xa. (2.1) 

with Xa homogeneous coordinates along X. 

We use a closed string version of the Berkovits model^ [11, 16] with a Euclidean world- 
sheet S. The fields are 

Z:T,^T, y G 0^'°(S) (g)r*T, and aGO°'^(S). 

^The standard Berkovits open twistor-string would be just as good for most of the considerations in this 
paper; although it is tied into split signature and gets some signs wrong, subtleties over choices of contours 
are avoided. For more on the version used here see [15]. 
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The action is 

S[Z, Y,a\ = J YiBZ^ + aZ^ Yj , (2.2) 
up to matter terms, and has the gauge freedom 

{Z,Y,a) ^ {e''Z,e-"Y,a-da), a G C. 

The gauging reduces the string theory to one in PT and the formaUsm allows one to use 
homogeneous coordinates on PT. 

Amplitudes are computed as an integral of worldsheet correlators of vertex operators 
on S over the moduli space of 'zero-modes': the space of classical solutions to the equations 
of motion. For gravity, the vertex operators correspond to deformations of the complex 
structure together with the cohomology class of deformations of the B-field. These are 
given by 5-closed (0, l)-forms F on the bosonic part of twistor space PT with values in 
the generalized tangent space T©T*PT. On T we represent these by 5-closed (0, l)-forms 
•= {f^t9i) of homogeneity (1,-1) satisfying djf^ = = Z^gj, defined modulo gauge 
transformations {aZ^ , djP). These conditions imply that {f^dj, gjdZ^) represents a section 
of T © r*Fir. The corresponding vertex operators take the form 



VF:=Vf + Vg:= J f{ZYYi+g{Z)idZ' . 



For n-particle tree-level amplitudes, we take S = CP^ and the amplitude reduces to 
an integral of a correlation function of n vertex operators 

oo „ 

M{1, ...,n) = J2 d/Xd(V>i . . . Vpjd , (2.3) 

over the space ^d,n of maps Z : CP^ PT of degree-ci and n marked points.^ This is 
the moduli space of zero-modes for the twistor-string, and has a mathematical definition 
as a supersymmetric generalization of Kontsevich's moduli space of stable maps [17]. To 
be more concrete, we can represent the maps by 

Z(.)=Y..i.t'Ur, dM. = ^;jQj!p;^nd'IU. (2.4) 

r=0 ^ ' ^ r=0 

where a a are homogeneous coordinates on CP"*^ and G T provide a set of coordinates on 
^d,o with redundancy GL(2,C) acting on the a and hence the Ur- The vertex operator 
Vp. = VF{Z{ai)) is inserted at the zth marked point (Tj G E, and the correlator naturally 

introduces a (l,0)-form at each marked point either from the Yj or the dZ^ , whereas the 
'wave-functions' {f^ ,gi) naturally restrict to give a (0, l)-form at each marked point. See 
[18] for further explanation. 



*The rules for taking the correlators are different at different degrees, hence the subscript d on the 
correlator. 
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The correlators are computed by performing Wick contractions of all the ys with Zs 
to give the propagator 

{Yia)jZ''{a'))d = (^)'^' ^ ' Oct = (a da) , (a a') = aoa[ - aia'^ , (2.5) 

where {aa') = eAsa^a'^ is the SL(2,C) invariant inner product on the world-sheet coor- 
dinates. When Y acts on a function of Z at degree d, it then differentiates before applying 
the contraction; Y acting on the vacuum gives zero so that all available Ys must be con- 
tracted, but this contraction can occur with any available Z. The ^ is an arbitrary point 
on the Riemann sphere and reflects the ambiguity in inverting the 9-operator on functions 
of weight d on S. The overall formula should end up being independent of the choice of ^. 

Unlike the Yang-Mills case, the degree d of the map is not directly related to the MHV 
degree of the amplitude (which essentially counts the number of negative helicity gravitons 
minus 2). The MHV degree of an amplitude counts the number of insertions of Vg minus 
2 (so the MHV amplitude has two V^s). Conformal supergravity amplitudes have been 
calculated from this formula in [19, 20]. 

2.2 Reduction to Einstein gravity 

We use a modified form of the Maldacena argument [10] in which the cosmological constant 
is not normalized, so that we can investigate the limit A — t- 0. The argument then shows 
that the classical S-matrix for conformal supergravity (i.e., the action evaluated perturba- 
tively on in- and out-wave functions) agrees with that for Einstein gravity, but multiplied 
by an additional factor of the cosmological constant A. 

Even for Einstein gravity without supersymmetry, we will need to use the superge- 
ometry of = 4 supertwistor space. ^ We first give the restriction required of the vertex 
operators for A/" = 4 supersymmetry and then for M = 0. 

To reduce to the Einstein case we must break conformal invar iance. This is done by 
introducing skew infinity twistors Ijj, I^^ with super-indices I,J = {a, a). The bosonic 
parts Ice/3, I"'^ satisfy 

where A is the cosmological constant and we shall set the fermionic part equal to zero. In 
terms of the spinor decomposition of a twistor Z°' = (A^, fi^') we have 




e' 



They have rank two when A = (i.e., the cosmological constant vanishes) and four oth- 
erwise. The fermionic parts of I^"^ can be non-zero and correspond to a gauging for the 
i?-symmetry of the supergravity [21, 22). 



■'At tree- level we can pull out the pure gravity parts of the amplitude, but their construction still 
relies on the fermionic integration built into the twistor-string formulae; this simply leads to some spinor 
contractions. 
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Geometrically /^"^ and //,/ respectively define a Poisson structure {, } of weight —2 and 
contact structure r of weight 2 by 

{hi, h2} := I^^dihidjh2 , T = IijZ^dZ^ , 

and we can use the Poisson structure to define Hamiltonian vector fields = I^'^djhdj 
which are homogeneous when h has weight 2. Lines in PT on which the contact form r 
vanishes correspond to points at infinity. This zero-set defines a surface ^ in space-time 
which is null when A = 0, space- like for A > 0, and time-like for A < 0. 

The Einstein vertex operators (V/j, VjJ correspond to Vf + Vg subject to the restriction 

{f,9i) = {l''djh,~hIijZ') 



so that 



[ I'^Yidjh, V^= / hAT. (2.7) 

The first of these gives the A/" = 4 multiplet containing the negative helicity graviton and 
the second that containing the positive helicity graviton. Sec [9] for more discussion. In 
order to reduce to standard non-supersymmetric Einstein gravity we must impose 

h = 62, and h = ■ (2-8) 

Thus, evaluated on vertex operators constructed from (2.7) with (2.8), (2.3) leads to the 
construction of Einstein gravity tree amplitudes. With this restriction, we see that there 
is now a correlation between degree of the maps and MHV degree, as fermionic variables 
only come with h: since there are 4d fermionic integrations in the path integral for the 
amplitude, there must be d insertions of h. 

At M = 4 there are spurious amplitudes that can be constructed from the Einstein 
wave functions corresponding to other conformal supergravity sectors. In [9], we argued 
that the Einstein gravity amplitudes could be isolated by imposing the correspondence 
between the degree of the map and the MHV degree: d = k + 1. This leads to the following 
starting point for tree-level scattering amplitudes in Einstein gravity from twistor-string 
theory: 



A1^(l, . . . , n) = lun i [ d;.,+i {v,, ■ ■ ■ V,„_,_,Vj^^_^_^ ■ ■ ■ I- 

= lim ^ / duk+i (y -dhi-'-Y ■ dh^-k-z K-k-2Tn-k-2 ■ ■ ■ Ktu , 

A /^^^^^ \ / fc+1 

= M 1 / d^^k^\ Cn,k+1 , (2.9) 

where Y ■ dhi = I^'^Yndj ihi and the last equation defines Cn,k+i as the worldsheet corre- 
lation function of the relevant vertex operators. 

In this paper we will focus on the A — >■ limit, in particular on the 0{A^) and 0(A) 
parts. However, according to the Maldacena argument, if the Berkovits-Witten formulation 
correctly gives conformal supergravity amplitudes, then (2.9) will do so for Einstein gravity 
for ah A. 
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3 Evaluating the Twistor-string Tree Formula 



We now turn to the evaluation of the formula (2.9) to test the claim that it will produce 
Einstein gravity amplitudes correctly. The most elementary consequence is that at A = 0, 
Cn,k+i must vanish because of the overall extra factor of A. Wc will also be interested in 
the coefficient of A since according to the version of the Maldacena argument, this will give 
the A = Einstein amplitudes which are now well known. 

3.1 Dependence on A 

Powers of A arise from the identity /"'^//j^ = AS". This gives a factor of A whenever there 
is a Wick contraction between a Y which is always be attached to an upstairs infinity 
twistor and a r which contains a downstairs infinity twistor as follows: 



= I '^^^ (e^ V ^ ^ ) 



where d' is the holomorphic exterior derivative with respect to a' . Thus, at 0{A^) we must 
exclude all Wick contractions with a r. At 0(A) we are allowed one such contraction and 
indeed we will see that the non-trivial contribution has precisely one such r-contraction. 

3.2 Naive computation of the 0{A^) contribution 

At 0{A^) we exclude any contraction of a 1^ with a r, so a 1^ can only contract with one of 
the wave functions. In the naive calculation we will allow all such contractions. However, 
we will see that this does not vanish even for the MHV amplitude. 

Working at A; = in (2.9), insert the standard momentum eigenstates 



hj 



[ - Pj) exp (isjiil^jXj]]) . (3.1) 



/C 

In the A — > limit we find 

{I^JYiidjhihj) = ^^^-^h,h,. 

After performing all world-sheet integrals, the naive contribution from the ith Yi to the 
A = part of the correlator is: 



E 



borrowing Hodges' notation [13]. Hodges shows that it follows from momentum conserva- 
tion that this is independent of ^ and indeed has an interpretation as an inverse soft factor 
for inserting a particle at the site i. This gives 

n-2 

Cn,l|o(AO) = n 



1=1 
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but this does not generically vanish. 



3.3 Connected trees and Feynman diagrams 

Rather than give up and deduce that the Berkovits-Witten twistor-string is incorrect, we 
instead use this as an indication that a different rule is needed for evaluating the twistor- 
string formula. The rule we will impose is that the Feynman diagrams for the contractions 
in the correlators appearing in (2.9) must be connected trees. This eliminates the loops 
and minimises the number of disconnected terms which are implicitly included in the naive 
calculation above. We will see that it leads to the Hodges [13] formula at MHV via the 
weighted Matrix- Tree theorem, an analogue of Kirchoff 's theorem for directed weighted 
graphs (see appendix A and [23-25]). It also sheds light on the possible twistor-string 
origin of Cachazo and Skinner's recent formula for tree-level N'^MHV amplitudes [5, 6]. As 
a bonus, it also provides an understanding of the origins for the tree formulae of [12] and 
[26] in twistor-string theory. 

Our argument represents terms contributing to the correlators of (2.9) as Feynman 
diagrams of the world-sheet CFT. It will be useful to represent such diagrams graphically, 
so we set out our notation here. There are two different kinds of vertices and two different 
kinds of propagator (or contraction) which can contribute to Cn,k+i- Vertices correspond 
to either Yi ■ dhi (white) or hjTj (black). Straight solid directed edges correspond to 
propagators arising from contractions directed outwards from a Yi to an Einstein wave 
function (either hj or hj); this utilizes the Poisson structure of the twistor space. Straight 
dashed edges correspond to a contraction directed outwards from a Yi to a Tj, and involve 
the contact structure. These both arise from standard contractions in the Berkovits-Witten 
twistor-string theory [11, 18]. See figure 1. We will only allow contributions from diagrams 

• kn O » O ^^^'^'^ = M wh) 

* i j 

O • (Y'^j) 

O Yrdhi i ^ 

i 

Figure 1. Building blocks for Feynman diagrams 

that are trees (i.e., no closed cycles) and connected after removing all black vertices. 

The diagrams acquire a factor of A*" if there are r {YiTj) contractions and we will just 
consider the cases where r = or r = 1.^ As we have said, the r = case should vanish 
and r = 1 with a single {Yt) contraction should give the Einstein amplitudes at A = 0; 
r > 1 will give contributions to the de Sitter amplitudes that are higher order in A. An 
example of a Feynman tree with r = 1 contributing to the five-point MHV amplitude is 
given in figure 2. 

We now investigate this formalism in detail for MHV amplitudes, and then generalize 
it to the N'^MHV case in a subsequent section. 

®There are also potential powers of A coming from the Poisson structure, but as we observe later, these 
do not contribute at 0(1) or 0(A). 



-9- 



5 



Figure 2. An example of a Feynman tree diagram contributing to 
3.4 The 0(A°) case 

For k = the diagrams will have 2 black and n — 2 white vertices. With A = there will 
be no {YiTj) contractions and the only edges in play are the directed ones from a vertex 
operator with a y to any twistor wave function. 

We begin by drawing a directed graph G which features all possible edges at once: 
each white vertex has n — 1 outgoing edges (see figure 3) to each of the other vertices, 
but just n — 3 incoming from the remaining white vertices. Each black vertex has just the 
incoming n — 2 edges from each of the white vertices. The (n — 2)(n — 1) edges of G can 
be labeled by a pair for the edge starting on white vertex i and ending on j ^ i. 

The Matrix- Tree theorem in its simplest form allows us to count the number of oriented 
spanning trees of G; later we will weight the counting by the contribution from all the 
contractions. The key device is the n x n Laplacian matrix of G constructed as follows. 

We first define the two modified incidence matrices, E and E for G with rows indexed 
by the vertices v oiG and columns by the edges e. The [v, e) entry of S is 1 if e is outgoing 
from V and otherwise, whereas E has +1 if e is outgoing and —1 if e is incoming and 
otherwise. For G, the components are given by: 



E. 



lifi = j 
otherwise 



E. 



1 if i = i 
-liii = k 
otherwise 



(3.2) 



The Laplacian matrix is then defined by 

/ n 



L = EE' 



1 -1 



V 



-1 o\ 



n 



-10 
-10 
-10 



(3.3) 



The Matrix- Tree theorem (see appendix A or [23-25] for some background) enables us to 
use L to count the Feynman trees entering into the 0{K^) part of the MHV correlator. 



Lemma 3.1 The number of non-trivial MHV Feynman tree diagrams at 0{hP) is 2 det Ln-i n 
2{n — l)"'~^, where Ln-i n is the matrix L with the n — andn^^ rows and columns deleted, 
and 1 subtracted from the diagonal. 
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n-l 




}n-l 



} n - 2 



'o: 



Figure 3. The graph G features all possible contractions within the 0(A°) class 



Proof: The MHV Feynman graphs are spanning trees of the graph G obeying our rules 
(so that, in particular, there is only one outgoing edge from each white vertex). By Euler's 
theorem for a tree the n — 2 edges connect n — 1 vertices, so one vertex is left out, and it 
must be one of the black ones: either n — 1 or n. Consider the class of graphs where n 
has no incoming edges: delete the n*'^ row and column from L and subtract one from the 
diagonal entries (as each white vertex no longer connects to vertex n) to obtain Ln- Now, 
let Ln be the Laplacian matrix of some new graph G' . The Feynman graphs in this class are 
the spanning trees of G' rooted at vertex n — 1; by the Matrix- Tree theorem (see appendix 
A), these are counted by det Ln-i n- By Cayley's theorem, this is equal to (n — The 
other class is counted in an equivalent fashion giving the same answer. □ 

To perform the sum of Feynman graphs in the correlator, we construct a weighted 
version £ of L which includes the propagators required by the Feynman rules for the 
contractions. We insert a matrix over the set of edges whose diagonal entries are the 
twistor-string propagators: 



P = diag((/.;.), 
so that we can define the weighted Laplacian as 



(3.4) 



C = EVE^ 



It. 



\ 



o\ 



T^^n-2 






(3.5) 



This differs from the matrix of (1.1) by a rank 2 error term, and from Hodges' 
matrix $ (the A; = form of (1.2)) by a conjugation and rank 2 error term: 



'^E = -T-'^^T+S, 



\ 



in-l 



T = diag ((C if) . 
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The matrix $ has co-rank 3 thanks to the relation 

n 

5^l>i,a/af = 0, (3.6) 

j=i 

which follows from momentum conservation [13]. Furthermore, the error term is such that 
once row and columns n — 1 and n are deleted, the matrices are conjugate. It is this high 
co-rank that will enable us to deduce that the sum of diagrams vanishes. 

Although Lemma 3.1 generalizes to tell us how to compute the relevant correlators by 
taking determinants of minors of £, we must divide up into sub-cases according to whether 
the remaining outgoing propagator from the white vertices lands on the black vertex n — 1 
or n. We first, consider the Feynman trees rooted at vertex n — 1; this means that there 
can be no directed edges to vertex n and we must delete from each diagonal entry. We 
therefore define 

Definition 3.1 Let £" be C with 4)\ deleted from the ith diagonal entry. We further define 
£" A; with rows and columns i, . . . ,k deleted. 

For diagrams not involving the vertex n, wc must also cross out the n^^^ row and column 
of C leaving the modified (n — 1) x (n — 1) matrix denoted C!^. The sum of correlators 
associated to diagrams in this class is then given by crossing out the row and column of 
corresponding to vertex n — 1 to yield and taking its determinant, thanks to 

the weighted version of the Matrix- Tree theorem (c.f., appendix A or [24], Chapter 36). A 
similar pattern follows for the other class of trees rooted at vertex n. This indicates that 
the A = portion of the conformal gravity MHV correlator is given by 

Cn,i|o(AO) = det/:;^„_i -)-deti:^Ii„. 

This sum can be re-interpreted thanks to some basic facts about the matrix £, and deter- 
minants. 

Lemma 3.2 We have 

detC:ij= ^ Rk<l>';, i,j,k = l,...,n, (3.7) 

k^...i,j 

where Rk = 

Proof: The matrix -C' '^^ only depends on the (fJ^ through its diagonal entries by the 
definition of C in (3.5). We can therefore expand its determinant out along the diagonal 

by 

det/::j^. = A+ Rk^^v 

where = so it suffices to show that ^4 = 0. By lemma 3.1 and the Matrix- Tree 

theorem, we know that each term in det£' ^ corresponds to the CFT correlator of a tree 
rooted at vertex j (omitting vertices . . . also) . This means that every term will contain (at 
least) one propagator to vertex j, by definition of a rooted tree graph. Thus there can be 
no term in det>C;;j ■ of order zero in <fH for all k ^ j; in other words, A = 0. B 
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Lemma 3.3 det „_i + det Clz\ „ = det Cn-i n- 
Proof: Lemma 3.2 means that, 

detCn-l= E I^l-\n= E ^^^n' 

i^n— l,n i^n— l,n 

where Ri and S'j arc the remaining factors. However, the matrices and >C"l} „ differ 

only by the presence of (l)\_x and along their diagonals, respectively, so it follows that 
Ri = Si- Therefore, we can write 

det + det £;::i„= E Rii<l>Li + ^i) = detCn-in, 

as required. □ 

Lemmas 3.2 and 3.3 therefore give 

det/:n-i„ = det(-(T-iir)n_i„ + ^„_i„) = det = 0, (3.8) 

since $ has co-rank three and the error term was only non- vanishing in the n — 1^* and n*^ 
columns. This confirms that the A = contribution to C„^i vanishes, as required by the 
Maldacena relation between conformal and Einstein gravity amplitudes. 



4 The MHV Amplitude From the 0(A) Contribution 

Hodges' formula for the MHV amplitude in A/" = 8 supergravity is [13]: 

M^ = J^I^^||x:^Ajdet'($°), (4.1) 

where det'($°) is a reduced determinant of the MHV {k = 0) matrix (1.2), and the delta- 
function expresses super-momentum conservation. Although other forms of the MHV am- 
plitude have been known for some time (e.g., the BGK or BCFW type formulas [27- 
29]), (4.1) manifests the full permutation symmetry of gravitational scattering amplitudes 
without a cumbersome summation, and makes no reference to any cyclic ordering of the 
gravitons. In this sense, Hodges' formula is the closest gravitational analogue to the Parke- 
Taylor formula for gluon scattering in Yang-Mills theory. It was proved using a M = 7 
version of BCFW recursion [30]. Here we will see how it arises from the twistor-string 
formula (2.9) at 0(A). 

There are two possible 0(A) contributions. The first arises from a Y-t contraction as 
we described in the previous section. The other arises from the 0(A) part of the Poisson 
structure and will potentially contribute to all possible contractions between a Y and a 
wavefunction. The second can be treated very similarly to the 0(A°) case and can be seen 
to have a vanishing contribution. 
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Figure 4. The graph showing all possible contractions subject to a given (Kir„) contraction 



We first of all ignore the rs so the contribution from the downstairs infinity twistor 
is irrelevant (the vanishing will be seen before any integration of a r is performed as 
in the previous section). The argument for the vanishing of the 0{A^) contribution to 
the correlator works for any choice of wave functions, so the calculation is essentially 
conformally invariant (we are working in a conformal gravity framework) and will work 
to show that the contribution of the connected tree diagrams with no Y-t contractions 
vanishes for any choice of simple upstairs infinity twistor. Thus the derivative of the 
correlator with respect to the infinity twistor vanishes in any direction for which the infinity 
twistor remains simple. But such directions span all possible variations in the infinity 
twistor: an infinity twistor corresponds a point in and the simple ones to lightcones 
in and any direction can be expressed as a sum of null directions. In particular, the 
derivative with respect to A must vanish and so the 0(A) term vanishes.^ 

Hence, the 0(A) contribution to the correlator is generated by those Feynman trees 
which have a single {YiTj) contraction. Without loss of generality, suppose this contraction 
occurs between vertex i and n. We will see that, remarkably, once the answer is added on 
to the case where Yi contracts with r„_i the result is independent of the choice of i. 

We proceed to draw a graph that contains all the remaining possible edges allowed 
by our Feynman rules (see figure 4). From this point, we can proceed almost identically to 
the 0(A'') case. The modified incidence matrices now read 



E. 



1 if j = k and j ^ i 
otherwise 



E. 



1 if j = k and j ^ i 
-lifj = l 
otherwise 



(4.2) 



and the weighted Laplacian matrix C = E^V{E^)'^ is the same as in (3.5) except that now 
the i*^ column has all zero entries. 

The contraction (YiTn) combined with (YiTn-i) produces a factor corresponding to the 
3-point MHV amplitude. We will discuss this factor in greater detail later, but we will take 



^An alternative argument using dual-twistor wave functions can be used to show directly that any 
choice of upstairs infinity twistor gives zero when Y-t contractions are ignored using precisely analogous 
arguments as the previous subsection but now with a Hodges like matrix in which the numerators are 
replaced by I^^WiaWjii and the co-rank 3 arises from an analogue of momentum conservation arising from 
the delta-functions in the Grassmannian approach of [6]. 
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its value^ (with Af = 4 supersymmetry) to be given by: 

) = ADaiDanDan-iy-- — vJT (4-3) 

[1 n)"'(ri — 1 zj^ 

Hence, we are left with the task of computing the factor arising from the remaining sum 
over trees in G*. 

There are n — 3 remaining l^s and these must connect the remaining n — 3 vertices 
into a tree so there will be one remaining outgoing contraction from this tree and this can 
contract with the vertex z, n — 1 or n. The arguments in the proof of lemma 3.1 show that 
there will be (n — 2)"~^ trees contributing, and that the sum of correlators is given by: 

det r L + det Clz'rn, + det C ^ „-i , (4-4) 

Here, ^ is the matrix C with the i**^, j*^, and k^^ rows and columns deleted, and the (p^, 
(pj terms subtracted from the remaining diagonal entries. A simple application of lemma 
3.2 tells us that 

detCrnZL = T.^^^n, detCl-_l';,, = ^Rj<H, det = Rj4i-i, 

j j 3 

and by lemma 3.3, this gives 

det „ + det Clzl , + det C ^ n-i = E ^^^^i + <^n-i + <l^n) = det A „-i „. 

3 

Noting that det £^ „_i „ = |. 

in-in, we then find that the correlator with 
a single contraction of the form {YiTn) gives a contribution: 

Q,ilo(A) = A ^^^),(^_^^), det^.„_i.. (4.5) 

Although wc still have to sum over the n — 2 different choices of i for the (YiTn) contraction 
(the interchange n -H- n — 1 is already taken care of in (4.3)), the defining properties of the 
reduced determinant det' of the matrix $ in fact show that all of these terms will be equal 
[13]. This is a consequence of the emergent permutation symmetry that arises in A/" = 8 
supergravity. 

Indeed, we can introduce a power of (n — 1 n)~^ into the three point MHV amplitude 
(4.3) by extending the supersymmetry to = 8 where now this emergent permutation 
symmetry will extend over n and n — 1 also. This gives the = 8 momentum-space 
amplitude: 

f d^^CV^d^" ~ detl>,„-i„ 

= volGL(2,C) nD-^-^^(^^-^-^(-^))(„„-l)2(n_i,)2(,„)2 

= <^^I^MX:A,A,jdet'(i°). (4.6) 



*We will see later that our nai've rules that we have used here give rise to a gauge-dependent error term. 
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Here, the super-momentum-conserving delta function is realized via Nair's lemma [31] 

^41 16 f^xX^ ^ J ^41 16^ (^iX^X^rx^^' + irjaXAO'' 



and det'($'^) is short-hand for: 

det'($") = 



//5.o^ _ det($i„_i„) 



{in — — 1 n)'^{n i)"^ ' 



with the GL(2,C) freedom fixed by setting Uj = Xj. 

Equation (4.6) is just a particular representation of Hodges' formula (4.1) in the form 
found by Cachazo and Skinner [5, 6], so our Feynman tree formalism correctly produces 
the MHV amplitude for gravity. 

4.1 The 3-point amplitude 

In arriving at (4.6), wc took for granted that our twistor-string theory correctly reproduced 
the three-point MHV amplitude with J\f = 4 supcrsymmctry: 

{in — l)^(n z)"^ 

where we have suppressed a momentum conserving delta-function. The Feynman diagram 
associated with the three point amplitude is depicted in figure 5, and the twistor-string 
expression for the correlator is given by: 



dm {Yi ■ dhi hn-lTn-l hnTn) . (4.8) 



However, the naive use of Wick contractions in the worldsheet CFT leads to a gauge 
dependent error term between (4.3) and (4.8) given by 

«'>«^> (4.9) 



This same error arises also if we directly use the three point formula given by Nair [18]. 

In [9] we instead used a geometric computation based on a technique utilized to con- 
struct the MHV amplitude in the Einstein case in [28]. This was based on the conformal 
gravity twistor action of [14] which closely models the twistor-string theory calculation. It 
is clear that the naive twistor-string rules that we have been using are not quite correct in 
this context given that they are leading to a gauge-dependent answer. 

In summary, if we assume that twistor-string theory produces the correct three-point 
MHV amplitude, then wc have a complete derivation of Hodges' MHV formula from twistor- 
string theory Feynman rules. However, there is a gap in our understanding of the three- 
point correlator in the world-sheet CFT and in particular the {Y r) contractions which are 
not yet properly understood. 
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Figure 5. One of the two graphs of the three-point amplitude {Yi ■ 9/iift-„T„/i„_iT„_i) 

5 The N^MHV Amplitude 

At N'^MHV, the Hodges MHV formula is generalised by the Cachazo-Skinner formula [5, 6]. 
Like the twistor-string formula, it is an integral over the moduli of rational curves of degree 
d = k — 1. It provides all tree-level N'^MHV amplitudes of A/" = 8 supergravity [6] and is 
given by: 



J vol GL(2, C) det^(^')det^(^') 11 (a. A.; . (5.1) 



Here, ^>'^ is the co-rank -|- 3 matrix (1.2), and is essentially its parity conjugate: 

The det'($'^) and det'(<I>'^) are modified reduced determinants, involving some additional 
Vandermonde determinant factors. In particular, det'(^'^) can be written 

~. _ det (l>f„_,_2...„) _ det (l-f^.fc.^...^^ 

~ n,-<fc,.>e{.,n-fc-2,...,n}(-?'^) " VM(z,n-fc-2,...,n)- ^^"^^ 

That (5.1) involves an integral over the moduli space of re-pointed stable maps of CP^ 
to twistor space is of course highly suggestive of a twistor-string origin. We will here not 
attempt to explain the whole formula, but just the factor det $f „_^_2...„ in terms of our 
tree graph arguments generalized from the MHV setting. The remaining terms should 
presumably also arise from the other ingredients in the twistor-string calculations but we 
have not yet successfully identified those terms. 

For the N'^MHV amplitude the twistor-string tree formula (2.9) 



liin \ I dfik+i Cn,k+i- 



A^O A 



is now an integral over the moduli of rational maps of degree d = k — 1 and there are now 
n — k — 2 vertex operators of the form Y ■ dh and k-\-2 oi the form hr. Again we will allow 
only maximally connected Feynman tree diagrams to contribute to the CFT correlator 
C-n,k+i] if we do not, then we will again see that we obtain incorrect non- vanishing answers 
at A = 0. The diagrams therefore now have k + 2 black vertices and n — k — 2 white 
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vertices. In particular there are just n — k — 2 YiS to be contracted with the hj, hj or tj and 
so the diagrams become correspondingly simpler with fewer edges. The number of edges 
corresponds to the size of the minors whose determinant we must take when we use the 
Matrix- Tree theorem. It turns out that the ranks of the matrices given in (1.1) and 
(1.2) encoding the twistor-string propagators are n — k — 3, decreasing as appropriate with 
k. 

Again at A = 0, where we must have vanishing contribution, there can be no {Yt) 
contractions and for the 0(A) contribution (which should give the Einstein amplitude) there 
should be just one. As before, we draw the master graph G of all possible Yi contractions 
with the hj and hj. Thus in G, black vertices have n — k — 2 incoming edges (one from 
each of the white vertices) and no outgoing edges, while white vertices will have n — k — 3 
incoming (one from each of the other white vertices) and n — 1 outgoing edges to all other 
vertices. The definition of the modified incidence matrices E, E and L = EE^^ for this 
new graph follows in precisely the same way as in the MHV setting. 

The incidence Laplace matrix of G is now: 



n 



1 _1 ... 



-1 



n- 1 • 
-1 • 



V 



-1 



-1 ••• 



with the last k + 2 columns composed entirely of zeros. For the propagators to form a 
connected tree, there can only be one outgoing propagator from the white vertices. If 
this lands on vertex n — k — 2, the other black vertices are isolated and the Matrix- Tree 
theorem gives the number of such Feynman trees rooted at the black vertex n — fc — 2 to 
be detL„_fc_2...„ = {n-k- l)"-'=-3 . 

To compute the actual correlator at 0(A°) rather than just the number of contributing 
trees, we define the weighted Laplacian C weighted by the appropriate propagators at degree 
d = k + 1 which gives 



-^L,_2 0---0\ 



in-k-2 



in—k—1 
^n-k-2 



^n-k-2 



••• 



••• 0/ 



k+2 



This C is related to by a conjugation by T = diag ((.^ i)^^"^^ and a rank k + 2 error 
term that disappears when the last n — k — 2 rows and columns are deleted. The natural 
generalizations of lemmas 3.2 and 3.3 then tell us that the sum of the Feynman trees 
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contributing to the correlator Cn,k+i is given by: 

detCn-k-2-n = det$^_fc_2...„ = 0, 

and this is vanishing because has co-rank k + 3. This gives the expected vanishing at 
A = 0. 

The Einstein gravity contribution should be given by the 0(A) contribution where 
there is just a single {YiTj) contraction, as in the MHV case. Instead of a 3-point MHV 

pre-factor, we expect to obtain something more complicated which encodes the data of the 
single contraction {YiTj), the rs and other Jacobian fators^. We will ignore such factors 
here and focus only from those that arise from contractions of the remaining Ys with the 
hs and ^s. 

As in the MHV amplitude, this amounts to trees rooted at one of vertex i or any of 
the k + 2 r- vertices. Yet again, counting these (n — A; — 2)""*^"^ trees corresponds to taking 
reduced determinants of the matrix C via the natural extensions of lemmas 3.1-3.3. Hence, 
we obtain the contribution det Hi 72— /c— 2---n 

= det ^in-k-2---n- coursc, we must still sum 
over all of the different ways of taking the initial (YiTj) contribution, but as in the MHV 
case, the properties of the matrix ensure that all of these contributions are equal up to 
the kind of Vandermonde factors that we are ignoring. 

6 Discussion 

Although our derivation is far from complete, we have nevertheless made significant progress 
in our check that the Berkovits-Witten twistor-string correctly reproduces the Einstein sub- 
sector of conformal gravity amplitudes. Now that we have the Cachazo-Skinner formula 
(5.1), our task has been to identify how terms in that formula can arise from twistor-string 
theory. We clearly have an explicit derivation of not only the moduli space integral and 
the external wave functions, but also the reduced determinant of in the formula. This 
arose by requiring that only connected tree diagrams should be taken to contribute to the 
CFT worldsheet correlators arising in the twistor-string amplitude formula; the reduced 
determinant arose by expressing the sum of the resulting Feynman diagrams as the deter- 
minant of a weighted Laplacian matrix for the graph of all possible contractions using the 
weighted Matrix- Tree theorem. Although this provides substantial, perhaps overwhelming, 
circumstantial evidence that the Berkovits-Witten twistor-string correctly gives the (tree- 
level) gravity sub-sector of conformal gravity, and hence circumstantial evidence that all 
of conformal gravity can be obtained in this way, we are still missing complete arguments 
even for a part (equation (4.9)) of the Hodges gravitational MHV formula. At this stage 
our argument does not reproduce any of det' $ nor the Vandermonde determinants in the 
N'^MHV Cachazo-Skinner formula. 

At this stage it seems likely that the Vandermonde determinants and det' ^> factors 
must all be present in the Berkovits-Witten twistor-string formula. However, a better 

®In the fc = 0, or d = 1 case, we were free to fix the GL(2,C) freedom so that ai = Ai; thereby setting 
the Jacobian factor equal to unity. For d > 1 this is no longer possible and a non-trivial factors can occur. 
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understanding of the worldsheet conformal field theory is required even for the three point 
amphtude. 

It is somewhat vexing that only tree diagrams are required in the worldsheet correlator. 
Prom the string theory perspective this is difficult to motivate as there are no parameters 
that can be used to suppress loop contributions. There may be some exponentiation argu- 
ment as in the Yang-Mills case [32]: there the twistor-string naively gives the determinant 
of the twistor 5-operator whereas the correct contribution to obtain Yang-Mills is the log 
of the determinant. However, the geometric interpretation of this restriction to trees is 
clear. The Y contractions are performing a perturbation expansion for the construction of 
holomorphic curves in a deformed twistor space given as solutions to 

daZ'^ix, a) = r^dph{Z) . (6.1) 

The Feynman trees build up the correct classical solution to this equation as appropriate to 
the Penrose non-linear graviton construction. The loop contributions would be appropriate 
to some quantisation of this construction, whereas all we are trying to obtain here are 
classical solutions in order to get gravity tree amplitudes (quantum ones according to 
string theory should only arise when higher genus curves are introduced). 

We have presented this work as a verification of the Berkovits-Witten twistor-string, 
but it is also the case that a key motivation is to elucidate the twistor-string underpinnings 
of the Cachazo-Skinncr formula. Although it is very suggestive of a twistor-string theory 
that directly describes A/" = 8 supcrgravity, the details arc far from clear. The ingredients 
described in this paper will be very likely to play a role also in such a theory, or at least in 
its A/" = 4 truncation, which should in turn embed in the conformal gravity twistor-string. 
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A Graph Theory 

In this appendix, we review some graph theory that we have used to sum the Feynman 
tree graphs. For a more complete introduction see e.g., [23-25]. 

Let G be a directed graph, with a set of vertices V and directed edges £. We will 
denote such a graph by G = (V, £). Denote the number of vertices as |V| = n and edges as 
\£\ = m. The incidence matrix of G is a nxm matrix E encoding the connective structure 
of the graph. We can denote an edge of G as (j, k) & £, meaning the directed edge which 
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is outgoing from vertex j G V and incoming to vertex k E V. In this notation, the entries 
of the incidence matrix are given by: 



Ei {j,k) 



' 1 if i = j 
-1 if ?; = /c . (A.l) 
otherwise 



The incidence matrix can be used to build the Laplacian matrix of the graph G, which 
is a n X n matrix given by: 



A = EE'^, a. 



■1 if z 7^ J and G S 

i{i=j , (A.2) 

otherwise 



where d°^^ is the out-degree of vertex i eV. 

The Feynman tree diagrams we are interested in counting throughout this paper are 
examples of particular sub-graphs of a directed graph, namely rooted spanning trees: 

Definition A.l ^4 spanning tree of G rooted at v G V is a spanning sub-graph T = {V,T) 
of G such that: (1.) T has no oriented cycles; (2.) vertex v has out-degree zero; and (3.) 
every other vertex i ^ v has out degree one. 

The main result in graph theory we use is the weighted generalisation of the Matrix- 
Tree theorem. The original matrix-tree theorem tells us how to count spanning trees of 
G rooted at a particular vertex in terms of determinants of the reduced Laplacian of G 
(c.f., [25], Chapter 9; [24], Chapter 36). Since it illustrates the enumerative nature of the 
arguments used in the body of the text, we provide a proof here: 

Theorem 1 (Matrix- Tree Theorem) Let k{G, v) denote the number of oriented span- 
ning trees of G rooted at v £ V, and Ay be the Laplacian matrix of G with the v^^ row and 
column deleted. Then k{G,v) = det A^. 

Proof: Without loss of generality, let f = n G V. Let the (n — 1) x n matrix F be 
the incidence matrix E oi G with its n**^ row deleted, so that A„ = FF"^ . Now, the 
Cauchy-Binet theorem (c.f., [25], Chapter 9) tells us that we can write: 



det A„ = ^ dei{Fs) det(Fj) = ^ dei{Fs) 



where S ranges over all subsets of £ of size n — 1 none of which arc outgoing from n G V, 
and Fs is the (n — 1) x (n— 1) matrix whose columns are those of F with indices in S. This 
means that each choice of S corresponds of a sub-graph T C G with n — 1 edges which 
have d°^^ = 0. We now quote the following lemma from graph theory: 

Lemma A.l (c.f., [25] Lemma 9.7) Let S be a set ofn—1 directed edges of G = (V,<?), 
none of which are outgoing from n G V. Then det Fs = ±1 if S forms a spanning tree of 
G rooted at n, and det Fg = Q otherwise. 
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This means that 



detA„= {±lf = K{G,n), 

spanning trees rooted at n 

as required. □ 

The weighted generahzations of the Matrix- Tree theorem include weighted edges (e.g., 
[24], Chapter 36) and allow us to introduce the propagator appropriate to each edge. This 
leads to the weighted Laplacians used throughout the paper. We define the weighted 
Laplacian matrix A for the directed graph G with weights for each edge given by: 

{-Wij ifiy^j and (i, j) G £ 
otherwise 

Then the weighted Matrix- Tree theorem reads: 

Theorem 2 (Weighted Matrix- Tree Theorem) Let T{G, v) he the set of all oriented 
spanning trees of G rooted at v ^ V and A be the weighted Laplacian of G. For any 
T G T{G,v), let £t denote its set of directed edges. Then 

det = ^ Yi ^ij 1 • 

TeT{G,v) \{i,j)&£T J 

When the weights are set to be twistor-string propagators, this is the primary result used 
in the text. 
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